UNCLASSIFIED 


Defense  Technical  Information  Center 
Compilation  Part  Notice 

ADPO 14041 

TITLE:  Doppler  Properties  of  Airborne  Clutter 

DISTRIBUTION:  Approved  for  public  release,  distribution  unlimited 
Availability:  Hard  copy  only. 


This  paper  is  part  of  the  following  report: 

TITLE:  Military  Application  of  Space-Time  Adaptive  Processing  [Les 
applications  militaires  du  traitement  adaptatif  espace-temps] 

To  order  the  complete  compilation  report,  use:  ADA415645 

The  component  part  is  provided  here  to  allow  users  access  to  individually  authored  sections 
of  proceedings,  annals,  symposia,  etc.  However,  the  component  should  be  considered  within 
the  context  of  the  overall  compilation  report  and  not  as  a stand-alone  technical  report. 

The  following  component  part  numbers  comprise  the  compilation  report: 

ADPO  14040  thru  ADPO  14047 


UNCLASSIFIED 


2-1 
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Tel  +49  228  9435  377;  Fax  +49  228  9435  618;  email  r.klemm@fgan.de 

Abstract 

Detection  of  slowly  moving  targets  by  air-  and  spacebome  MTI  radar  is  heavily  degraded  by  the  motion  induced 
Doppler  spread  of  clutter  returns.  Space-time  adaptive  processing  (STAP)  can  achieve  optimum  clutter  rejection 
via  implicit  platform  motion  compensation.  In  this  report  the  fundamentals  and  properties  of  STAP  applied  to 
air-  and  spacebome  MTI  radar  are  summarised.  The  effect  of  platform  motion  on  the  characteristics  of  airborne 
clutter  is  discussed.  The  performance  of  the  optimum  space-time  processor  is  shown.  Comparison  with  spatial 
or  temporal  only  techniques  illustrate  the  importance  of  space-time  processing. 

1 Introduction 

The  main  application  of  space-time  adaptive  processing  (STAP)  is  the  suppression  of  clutter  received  by  a 
moving  radar.  In  this  case  we  talk  about  space-slow  time  (pulse-to-pulse)  processing.  The  radar  platform 
motion  causes  clutter  returns  to  be  Doppler  shifted.  The  Doppler  shift  is  proportional  to  the  platform  velocity 
and  the  angle  of  arrival.  The  total  of  all  clutter  arrivals  sums  up  in  a Doppler  broadband  clutter  echo.  Targets 
whose  Doppler  frequencies  fall  into  the  clutter  Doppler  bandwidth  may  be  difficult  to  detect.  It  should  be 
noted  that  most  moving  targets  are  ’’slow”  or  low  Doppler  targets,  either  because  they  are  slow  or  exhibit  a low 
Doppler  due  to  their  motion  direction.  For  a radar  on  a satellite  (10  km/s)  all  targets  near  the  ground  (jet  aircraft 
300  m/s)  are  slow  targets. 

Other  important  techniques  operate  in  the  space-fast  time  (range  equivalent)  domain.  They  are  used  to 
either  suppress  jammers  in  broadband  array  radar  or  for  mitigating  terrain  scattered  jamming. 

Space-time  processing  needs  space-time  data  to  operate  on.  Space-time  data  are  obtained  from  a radar 
which  has  a phased  array  antenna  with  multiple  outputs  (spatial  dimension),  and  transmits  coherent  pulse  trains 
(temporal  dimension). 

1.1  The  principle  of  adaptive  clutter  suppression 

1.1.1  Practical  application:  adaptive  clutter  filter  fpr  surveillance  radar 

The  following  experiment  conducted  in  the  early  1970’s  (BDhring  & KLEMM  [6])  is  decribed  briefly  to  il- 
lustrate the  principle  of  adaptive  clutter  filtering.  This  filter  was  designed  to  suppress  weather  clutter  with 
unknown  centre  Doppler  frequency  and  bandwidth.  This  adaptive  MTI  systems  was  operated  with  a conven- 
tional grouhdbased  surveillance“=radar  with  rotating  reflector  antenna.  Figure  1 shows  the  block  diagram  of- 
an  adaptive  FIR  filter  based  on  the  prediction  error  filter  principle.  The  filter  coefficients  were  estimated  in 
real-time  and  the  echo  data  were  clutter  filtered  during  the  following  revolution  of  the  antenna. 

In  Figure  2 the  filtering  effect  is  demonstrated  using  simulated  clutter.  The  picture  shows  a photograph  of  a 
PPI  (pla  position  indicator)  screen.  The  zero  Doppler  filter  has  been  switched  off  so  that  a lot  of  ground  clutter 
can  be  noticed.  The  spokes  are  the  simulated  clutter  (for  simplicity  this  clutter  was  simulated  independent 
of  range).  The  available  radar  operated  in  L-band  which  is  quite  insensitive  to  weather.  Therefore,  to  have 
reproduceable  clutter  conditions  a hardwired  simulator  was  developed.  The  clutter  filter  was  adapted  based  on 
clutter  data  in  the  window  on  the  right.  As  can  be  seen  the  clutter  has  been  removed  and  a (true)  target  is  visible. 

Figure  3 and  4 show  suppression  of  a real  weather  cloud  before  and  after  filtering.  Except  for  a few  false 
alarms  the  weather  clutter  has  been  removed.  The  air  targets  (big  spots  on  the  left)  did  not  do  us  the  favour  to 
enter  our  measurement  window,  they  just  bypass  it1 . 

'No  responsible  pilot  will  enter  a thunder  storm  deliberately. 


Paper  presented  at  the  RTO  SET  Lecture  Series  on  "Military  Application  of  Space-Time  Adaptive  Processing  ", 
held  in  Istanbul,  Turkey,  16-17  September  2002;  Wachtberg,  Germany,  19-20  September  2002; 
Moscow,  Russia,  23-24  September  2002,  and  published  in  RTO-EN-027. 
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1.2  The  History  of  STAP 

In  the  following  two  sections  a brief  review  of  the  history  of  research  and  experimental  work  on  STAP  is  given. 

1.2.1  Theoretical  work  (Table  1) 

Table  1 includes  a number  of  publications  which  represent  milestones  in  the  evolution  of  STAP.  This  selection 
of  papers  is  certainly  not  exhaustive  but  to  our  opinion  representative.  In  some  cases  merely  the  earliest  papers 
rather  than  the  most  significant  papers  are  quoted.  'In  the  meantime  special  journal  issues  (#  20,21)  and  tutorials 
(#  12,19,20)  on  STAP  are  available. 

The  original  idea  of  reducing  the  clutter  spectrum  via  motion  compensation  originates  from  ANDERSON  (# 

1) .  This  paper  was  written  before  digital  technology  was  available.  The  DPCA  technique  Skolnik  (#  2)  was 
also  originally  invented  for  use  with  for  implementation  in  microwave  technology.  This  technique  compensates 
physically  or  electronically  for  the  platform  motion  in  a non-adaptive  fashion.  The  paper  by  Brennan  & 
Reed  (#  3)  is  the  basis  for  all  future  space-time  processing.  It  deals  with  interference  rejection  in  broadband 
arrays.  This  principle  was  analysed  in  some  detail  by  COMPTON  (#  7).  The  first  paper  on  STAP  was  written 
by  BRENNAN  et  al.  (#  4).  Here  the  optimum  (maximum  likelihood)  processor  applied  to  clutter  rejection  for 
moving  radar  has  been  described  and  analysed.  Starting  from  (#  4)  KLEMM  discovered  that  the  size  of  the 
clutter  subspace  is  about  N + M under  certain  conditions  which  gave  rise  to  manyfold  research  on  subspace 
processor  architectures  (#  8, 9,  12,  13,  18, 20).  Later  on  an  extension  of  this  result  became  known  as  ’’Brennan’s 
Rule”  (#  10).  Problems  of  real-time  implementation  of  STAP  processors  have  been  discussed  by  Farina  et 
al.  (#  11).  Ward  (#  15)  presented  angle  and  Doppler  estimation  errors  for  STAP  radar.  ENDER  (#  17) 
demonstrated  the  detection  and  re-positioning  of  moving  targets  in  SAR  images  obtained  with  the  multi-channel 
SAR  AER II.  DOHERTY  et  al.  and  JOUNY  et  al.  (#  16)  used  STAP  techniques  to  mitigate  the  effect  of  terrain 
scattered  jamming.  STAP  in  conjunction  with  bistatic  radar  has  been  discussed  by  Klemm  (#  23).  A recent 
trend  of  moving  reconnaissance  function  from  airborne  to  space-borne  platforms  can  be  noticed  (CovaULT 
(#  24)).  Besides  the  STAP  research  activities  in  the  USA  and  Europe  some  considerable  interest  of  Chinese 
scientists  in  STAP  can  be  noticed  (#  6). 

1.2.-2  Experiments  and  Systems  (Table  2) 

In  Table  2 experimental  and  operational  STAP  systems  are  listed.  The  first  (non-adaptive)  DPCA  experiment 
involving  an  array  antenna  has  been  carried  out  by  Tsandoulas  (#  1).  The  NRL  (#  3)  and  MCARM  experi- 
ments (#  4)  use  linear  sidelooking  arrays.  Both  programs  lead  to  many  detailed  investigations  on  possibilities 
of  airborne  clutter  rejection.  A large  number  of  publications  originate  frong these  programs,  dealing  with  var- 
ious research  topics  such  as  clutter  homogeneity,  sidelobe  Doppler  clutter,  knowledge-based  MTI  processing, 
subspace  techniques-(e.g.,  S — A),  and  bistatic  operation  (#  11).  The  Mountaintop  program  (#  5)  was  started  in 
1 990  to  study  advanced  processing  techniques  and  technologies  to  support  the  requirements  of  future  airborne 
early  warning  radar  platforms.  In  particular  the  effect  of  terrain  scattered  jamming  has  been  studied.  In  Europe 
the  AER  II  program  (#  7,  Germany),  the  DO-SAR  experiment  (#  8,  Germany)  and  the  DERA  experiment  (#  9, 
UK)  have  been  conducted. 

There  are  three  operational  systems  with  space-time  ground  clutter  rejection  capability  (Joint  STARS  (# 

2) ,  AN/APG-76  (#  6),  and  the  AN/APY-6).  The  first  one  has  a 3-aperture  sidelooking  array  antenna  and  has 
been  flown  in  the  Gulf  War.  The  AN/APG-76  is  a forward  looking  nose  radar  and  the  AN/APY-6  has  both 
sidelooking  and  forward  looking  capability.  From  the  available  literature  it  is  not  obvious  whether  these  system 
are  based  On  adaptive  algorithms  (STAP)  or  used  some  non-adaptive  DPCA-like  techniques. 

1.2.3  Historical  note 

The  non-adaptive  DPCA  technique  has  been  described  already  in  SKOLNIK  ed.  [39].  Research  on  space -time 
adaptive  processing  started  with  the  paper  by  BRENNAN  et  al.  [4]  in  1976  on  space-time  MTI  (moving  target 
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indication)  processing  for  airborne  radar.  Very  little  has  been  published  in  the  following  years.  In  1983  the 
author  introduced  the  concept  of  eigenanalysis  of  the  space-time  clutter  covariance  matrix  which  opened  the 
horizon  towards  subspace  techniques  for  real-time  applications  [21].  A number  of  follow-up  papers  by  the 
author  have  been  concerned  with  various  suboptimum  approaches  based  on  order  reducing  transforms  of  the 
signal  subspace.  These  papers  have  been  summarized  in  a book  on  STAP  (Klemm  [26]).  Since  1990  research 
activities  on  STAP  increased  tremendously,  particularly  in  the  USA  (e.g.  Ward  [46],  Wang  & Cai  [43]),  in 
China  (e.g.  Wang  & Bao  [44]),  in  the  UK  (Richardson  & Hayward  [35]),  and  in  Italy  Farina  et  al. 
[15]). 


# 

year 

subject 

authors 

1 

1958 

First  paper  on  motion  compensated  MTI 

Anderson  [1] 

2 

1970 

DPCA 

Skolnik  [39] 

3 

1973 

Theory  of  adaptive  radar 

Brennan  & Reed  [3] 

4 

1976 

First  paper  on  STAP 

Brennan  et  al.  [4] 

5 

1983 

Dimension  of  clutter  subspace 

Klemm  [21] 

6 

« 1987 

First  Chinese  papers 

e.g.  Wang  & Bao  [44] 

7 

1988 

Broadband  jammer  cancellation 

Compton  [9,  8] 

8 

1990 

Space-time  FIR  filter 

Klemm  & Ender  [22] 

9 

1992 

Spatial  transform  techniques 

Klemm  [23] 

10 

1992 

Brennan’s  Rule 

Brennan  & Staudaher  [5] 

11 

1992 

Real-time  implementation  of  STAP 

Farina  et  al.  [15] 

12 

1994 

Report  on  STAP 

Ward  [46] 

13 

1994 

Beam/Doppler  space  processing 

Wang  & Cai  [43] 

14 

1995 

STAP  for  forward  looking  arrays 

Richardson  & Hayward  [35] 
Klemm  [24] 

15 

1995 

Angle/velocity  estimation  with  STAP  radar 

Ward  [47] 

16 

1995 

Mitigation  of  terrain  scattered  jamming 

Doherty  [11],  Jouny  et  al.  [20] 

17 

1996 

STAP  for  SAR 

Ender  [13] 

18 

1996 

£ — A STAP 

Wangetal.  [45] 

19 

1996 

Effect  of  platform  maneuvers 

Richardson  et  al.  [36] 

20 

1998 

IEE  Colloquium  on  STAP 

Klemm  (chair)  [25] 

21 

1998 

Textbook  on  STAP 

Klemm  [26] 

22 

1998 

Effect  of  range  walk 

Kreyenkamp  [33] 

23 

1999 

IEE  ECEJ  special  issue  on  STAP 

Klemm  (ed)  [28] 

24 

1999 

IEEE  Trans.  AES:  Special  issue  on  STAP 
and  Adaptive  Arrays 

Melvin  (ed.)[32] 

25 

1999 

STAP  with  bistatic  radar 

Klemm  [27] 

26 

1999 

STAP  for  future  observation  satellites 

Covault  [10] 

Table  1.  Some  milestones  in  STAP  research 


echo  data 


Figure  1 : Block  diagram  of  adaptive  temporal  clutter  filter 


# 

year 

system 

authors 

1 

1973 

DPCA  with  array  antenna  (flight  experiment  ) 

Tsandoulas  [42] 

2 

1991 

Joint  STARS  antenna 

Shnitkin  [38] 

3 

1992 

NRL  experiment 

Lee  & Staudaher  [34] 

4 

1994 

MCARM  experiment 

Babu  [2] 

5 

1994 

Mountaintop 

Titi  [41] 

6 

1996 

AN/APG-76  sidelooking  airborne  array  radar 

Tobin  [40] 

7 

1996 

AER1I 

Ender  [12] 

8 

1996 

Domier’s  DO-SAR 

Hippier  & Fritsch  [19] 

9 

1996 

DERA  STAP  experiment 

Coe  et  al.  [7] 

10 

1998 

AN/APY-6  airborne  array  radar 

Gross  & Holt  [16] 

11 

1999 

STAP  with  bistatic  radar 

Sanyal  et  al.  [37] 

Table  2.  Existing  STAP  systems 


Figure  3:  Weather  clutter  before  adaption 
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2 Principle  of  air-  and  spaceborne  MTI  radar 

In  this  section  some  important  features  of  airborne  clutter  echoes  are  briefly  discussed.  The  efficiency  of  space- 
time  clutter  suppression  depends  significantly  on  these  properties.  For  definition  of  the  geometry  see  Figure  5. 
It  shows  two  important  cases  (sidelooking  and  forward  looking  arrays). 

2.1  Effect  of  Platform  Velocity 

2.1.1  Models  of  clutter  and  target 

The  results  presented  in  this  paper  have  been  calculated  on  the  basis  of  simple  models  for  target  and  interference. 
For  the  sake  of  brevity  we  give  here  only  the  models  as  used  in  the  evaluation.  For  more  details  the  reader  is 
referred  to  [26,  chapter  2], 


Clutter 

The  NM  x NM  clutter  covariance  matrix  has  the  form 


^ Qu 

Ql2  • • • 

Qim 

\ 

Q21 

Q22  • • • 

Q2  M 

Q ~ 

\ Qa/1 

QjW2  • • • 

Qmm 

/ 

where  the  indices  of  the  submatrices  m,p  denote  time  (pulse  repetition  intervals)  while  the  indices  i,  k run 
inside  the  submatrices  and  denote  space  (sensors).  N is  the  number  of  sensors  and  M the  number  of  coherently 
processed  echoes.  The  elements  of  Q are  integrals  over  a full  range  circle 


ll\n  — Pc  f PrnpPi—k  (2) 

J<p=  0 

xD2(ip)L2{ip)G{if,m)Gr{ip,p) 

X^m-piV,  vP)$\-k(<p)dp  + Pn 
m,p—l...M  i,  k = 1 . . . N 


where  Pc  is  the  clutter  power  at  the  single  element  at  a certain  instant  of  time  and  jfy  the  receiver  noise  power. 
The  other  symbols  denote  as  follows:  ip  azimuth,  s.  Figure  5;  A complex  clutter  amplitude;  £>(</?)  sensor 
directivity  pattern;  L(<p)  clutter  reflectivity;  spatial  (sensor-to-sensor)  correlation  which  is  not  considered 
here;  pmp  temporal  (echo-to-echo)  correlation;  G (<p,  m)  transmit  directivity  pattern.  The  temporal  and  spatial 
phase  terms  are  as  follows 


(v7)  = exp[/  —2 vpmT  cos  p cos  9] 

A 

2n 

= expb’  -y  {xi  cos  ip  + yi  sin  <p)  (3) 

x cos  9 - Zi  sin  9] 


The  indices  (,n  of  the  NM  x NM  covariance  matrix  are  related  to  the  sensor  indices  i,k  and  echo  indices 
m,p  through 


l — (rn  — l)iV  + i in  — l . 

. . M;  i — 1 . . 

..N 

(4) 

7 1 = (p—  1 )N  + k p = 1 . . 

t-H 

II 

..N 

(5) 

2-9 


Noise 


Receiver  noise  is  assumed  to  be  uncorrelated  in  space 

and 

time 

PiPm-Wp}  — ^ 

Pn 

0 

: m=p 

: m^p 

(6) 

E{THTirk}  = | 

Pn 

0 

: i = k 
: z k 

(7) 

where  Pn  denotes  the  white  noise  power. 


Target 

A target  at  azimuthal  position  <pt  moving  at  a radial  velocity  tyad  produces  the  following  space-time  signal  at 
the  array  output 

27r 

Smi(<Pt)  = Aexp[j—  (8) 

x (2urarf?n21  + (xi  cos  tpt  + yi  sin  ipt)  cos  9 
— Zjsin#)] 

m = 1 . . . M i — 1 ...  N 

2.1.2  Thelsodops 

Surfaces  of  constant  Doppler  are  given  by  cones.  Intersections  of  such  cones  with  the  planar  ground  results  in 
a set  of  hyperbolas  as  shown  in  Figure  6 for  horizontal  flight  and  flat  earth.  The  radar  platform  moves  from  left 
to  right.  Maximum  Doppler  is  encountered  at  0°  (positive  Doppler)  and  180°  (negative  Doppler),  zero  Doppler 
at  90°  270°. 

2.1.3  Impact  of  Array  Geometry 

Only  one  quadrant  of  the  isodop  field  is  shown  in  Figure  7 (thin  curves). 

The  axis  of  a linear  array  in  sidelooking  orientation  coincides  with  the  flight  path.  Therefore,  curves  of 
constant  look  direction  are  again  hyperbolas  on  the  ground  (fat  curves).  The  beam  traces  of  a linear  sidelooking 
array  coincide  with  the  isodops.  This  means  that  the  clutter  Doppler  is  range  independent.  This  important 
property  has  implications  on  clutter  rejection. 

The  axis  of  a forward  looking  linear  array  is  perpendicular  to  the  flight  axis.  Therefore,  the  set  of  beam 
traces  are  rotated  by  90°,  s.  Figure  8.  Now  one  can  notice  that  beam  traces  and  isodops  cross  each  other 
frequently.  For  a forward  looking  array  the  Doppler  frequency  of  clutter  echoes  depends"  obviously  on  range. 
The  range  dependency  occurs  especially  at  short  range,  that  is,  where  the  range  is  of  the  order  of  magnitude  of 
the  platform  height  above  ground. 

2.1.4  Azimuth-Doppler  Clutter  Trajectories 

Figure  9 shows  the  trajectories  of  clutter  spectra  in  the  azimuth2  (abszissa)  - Doppler  (ordinate)  plane.  The  four 
plots  have  been  calculated  for  different  crab  angles  if).  For  ^ — CP  (sidelooking  array)  all  the  clutter  power  is 
located  on  the  diagonal  of  the  plot  while  for  increasing  crab  angle  one  obtains  ellipses.  Finally,  ^ = 9(?  means 
a forward  looking  array,  with  the  clutter  power  located  on  a circular  trajectory. 

Clutter  echoes  depend  in  general  on  range.  First  of  all  the  backscattered  clutter  power  decreases  with  range 
according  to  the  radar  equation:  Pc  oc 

2cosy 
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Doppler-range  depence 

The  range  dependence  of  the  clutter  Doppler  frequency  has  already  been  addressed  in  the  context  of  isodops. 
It  follows  from  the  above  considerations  that  the  clutter  Doppler  is  constant  with  range  for  a sidelooking  array. 
Therefore,  all  four  curves  calculated  for  different  ranges  coincide  (Figure  9a).  From  Figure  9d  it  is  obvious  that 
for  a certain  look  direction  the  Doppler  frequency  increases  with  range  in  case  of  a forward  looking  array3. 

2.2  Comparison  of  Spatial,  Temporal  and  Space-Time  Processing 

The  principle  of  space-time  adaptive  processing  for  clutter  rejection  in  moving  radars  is  illustrated  in  Figure 
1 1 . A sidelooking  sensor  configuration  was  assumed.  The  clutter  spectrum  extends  along  the  diagonal  of  the 
cosp-up  plot.  Notice  the  modulation  by  the  transmit  beam. 

Conventional  temporal  processing  means  that  the  projection  of  the  clutter  spectrum  onto  the  up  axis  is 
cancelled  via  an  inverse  filter.  Such  filter  is  depicted  in  the  back  of  the  plot.  As  can  be  seen  the  clutter  notch  is 
determined  by  the  projected  clutter  mainlobe  which  is  a Doppler  response  of  the  transmit  beam.  Slow  targets 
are  attenuated. 

Spatial  processing  as  being  used  for  jammer  nulling  requires  that  the  clutter  spectrum  is  projected  onto  the 
cos  p axis.  Applying  an  inverse  spatial  clutter  filter,  however,  forms  a broad  stop  band  in  the  look  direction  so 
that  the  radar  becomes  blind.  Both  fast  and  slow  targets  fall  into  the  clutter  notch. 

Space-time  processing  exploits  the  fact  that  the  clutter  spectrum  is  basically  a narrow  ridge.  A space- 
time  clutter  filter,  therefore,  has  a two-dimensional  narrow  clutter  notch  so  that  even  slow  targets  fall  into  the 
passband. 


3This  follows  from  the  depression  angle  term  cos  9 in  (3). 
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Figure  9:  /r  — cos  <p  clutter  trajectories  for  linear  ar- 
rays: a.  ip  — 0°;  b.  ip  — 30°;  c.  ip  = 60°;  d.  ip  = 90°; 
from  inside  to  outside:  R/H  = 1.5;  2;  2.5;  3 


Figure  11:  Principle  of  space-time  clutter  filtering 
(sidelooking  array  antenna) 


Figure  10:  Range  dependence  of  the  clutter  Doppler 
frequency  for  forward  looking  array:  + p — 9(P  (look 
direction=flight  direction);  x /?  = 60° ; * fi  = 30° 
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3 Characteristics  of  Air-  and  Spaceborne  Clutter 

3.1  The  Space-time  Covariance  Matrix 

The  space-time  covariance  matrix  was  defined  in  eqs.  1,  2.  Figs.  12, 13  show  the  modulus  of  typical  space-time 
covariance  matrices.  As  can  be  seen  from  Fig.  12  the  spatial  submatrices  are  unity  matrices  with  the  diagonal 
shifted  with  the  temporal  indices  M,  p.  In  case  of  pure  spatial  processing  we  would  deal  with  a N x N 
unity  matrix  only.  In  this  case  no  substantial  gain  in  clutter  rejection  can  be  achieved.  Through  spacetime 
processing  we  obtain  the  other  correlation  ridges  in  the  matrix  which  provide  the  correlation  required  for  clutter 
cancellation.  If  we  use  directive  sensors  and  a directive  transmit  array  some  additional  correlation  comes  up  as 
can  be  seen  in  Fig.  13.  Fig.  14  shows  a space-time  covariance  matrix  for  noise  jamming.  As  can  be  seen  there 
is  no  temporal  correlation. 


3.2  Clutter  Spectra 
3.2.1  Eigenspectra 

The  concept  of  eigenanalysis  of  space-time  clutter  covariance  matrices  was  introduced  by  the  author  [21].  The 
eigenspectrum  (rank  ordered  sequence  of  eigenvalues)  shows  how  large  the  clutter  subspace  is. 

The  elements  of  the  space-time  covariance  matrix  are  calculated  as 

(Iln  ~ E{cimckp}  "b  Pu^ikrnp  (9) 

where  Cjm  is  the  space-time  clutter  signal.  The  spatial  (sensor)  indices  i,  k and  the  temporal  (echo  pulse)  indices 
m.p  are  related  with  the  matrix  indices  through 

l = (m-l)N  + i m = l...M;  i = l...N  (10) 

n = (p  — 1)N  + k p=l...  M\  k = l...N  (11) 

Pu,  is  the  receiver  noise  power  and  the  Kronecker  symbol.  It  was  found  in  [21]  that  for  a sidelooking 
equidistant  array  and  the  PRF  chosen  so  that  DPCA  conditions  are  fulfilled  (see  section  4.1.1)  the  number  of 
eigenvalues  is 

Ne  — N + M-  1. 

This  figure  determines  the  minimum  size  of  the  number  of  degrees  of  freedom  of  the  space-time  clutter  filter.4 


3.2.2  Power  Spectra 

Based  on  the  space-time  clutter-taoise  covariance  matrix  azimuth-Doppler  spectra  can  be  generated,  either  by 
2D  Fourier  transform  of  the  covariance  matrix,  or  by  use  of  one  of  the  well-known  high  resolution  power 
estimators. 

Let  us  consider  a covariance  matrix  of  the  form 


R = £{xx*}  = S -+-  N (12) 

where  N is  the  noise  component  and  S includes  all  kind  of  signal  or  interference.  Then  the  output  of  a signal 
matched  filter  is  simply 

Vsm{Q)=**  s(0)  (13) 


and  the  normalised  power  output  is 


Psm(&)  = 


s*(0)Rs(0) 

s*(0)s(0) 


(14) 


4 Actually  in  [21]  it  reads  Ne  — N + M.  The  correct  number  is  Nc  = N + M . 1. 
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s(0)  is  a steering  vector  which  seeks  for  signal  components  s(Q)  in  R.  P(&)  attempts  to  become  maximum 
wherever  the  steering  vector  5(0)  coincides  with  a signal  vector  s(6j)  in  R.  For  sinusoidal  signals  the  signal 
matched  filter  becomes  the  2D  Fourier  transform. 

Fig.  18  shows  a 2D  Fourier  clutter  spectrum.  One  recognises  the  main  beam  response  and  the  sidelobe  re- 
sponse along  the  diagonal.  In  addition  there  are  spurious  Doppler  and  azimuthal  sidelobes  which  are  responses 
of  the  spatial  and  temporal  FT  to  the  main  beam  clutter.  Notice  that  only  the  clutter  along  the  diagonal  is  phys- 
ical clutter,  the  sidelobes  along  the  Doppler  and  azimuth  axes  are  artifacts.  Similar  relations  can  be  observed  if 
Hamming  weighting  is  applied  in  space  and  time  (Fig.  19). 

The  minimum  variance  estimator  has  proven  to  be  the  most  useful  because  its  response  is  closest  to  the 
clutter  contained  in  the  covariance  matrix.  In  analogy  with  (19)  the  minimum  variance  estimator  becomes 

wmv  = 7R_1s  (15) 

with  7 = (s’R^s)-1.  The  output  power  is 

Pmv(&)  = (s*(0)R_1s(0))_1  (16) 

Figures  20  and  21  show  typical  MV  clutter  spectra  (F  is  the  normalised  Doppler  frequency)  for  sideways  and 
forward  looking  linear  arrays.  For  the  sidelooking  array  the  clutter  power  is  distributed  on  the  diagonal5  of  the 
plot  while  in  case  of  a forward  looking  configuration  the  clutter  is  distributed  on  a semi-circle.  Notice  that  in 
this  example  the  individual  sensors  have  sensor  patterns  which  blank  the  rear  semi-plane  out.  Otherwise  the 
clutter  power  of  the  forward  looking  array  would  be  distributed  on  a full  circle. 


5The  wavy  shape  of  the  trajectory  is  due  to  the  fact  that  we  used  tp  instead  of  cos  <p  for  the  abszissa. 


Figure  12:  Modulus  of  the  clutter  covariance  matrix 
vs  horizontal  and  vertical  index  (sidelooking  linear  ar- 
ray, N=12,  M=5,  omnidirectional  sensors  and  trans- 
mission) 


Figure  13:  Modulus  of  the  clutter  covariance  matrix  vs 
horizontal  and  vertical  index  (sidelooking  linear  array, 
N=12,  M=5,  directive  sensors  and  transmission) 
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Figure  14:  Jammer  covariance  matrix  (sidelooking  ar- 
ray, N=12.  M=5,  1 jammer,  absolute  values) 


Figure  15:  Eigenspectra  of  a spatial  covariance  matrix. 
N — 72,  M = 1;  o omnidirectional  sensors  and  trans- 
mission; * directive  sensors?-  x directive  sensors  and 
transmission 
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Figure  16:  Eigenspectra  of  a space-time  covariance 

matrix  (sidclooking  array).  N — 24,  M = 3;  o omni-  figure  18:  Fourier  clutter  spectrum  (sidelooking  array, 
directional  sensors  and  transmission;  * directive  sen-  <PL  = 9®  ) 
sors;  x directive  sensors  and  transmission 


Figure  17:  Eigenspectra  of  a space-time  covariance  Figure  ]9;  Fourjer  dutter  spectrum  with  spatial 
matrix  (forward  looking  array).  N — 24,  M = 3;  o ancj  temporal  Hamming  weighting  (sidclooking  array, 
omnidirectional  sensors  and  transmission;  * directive  _ o°) 
sensors;  x directive  sensors  and  transmission 
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4 The  Optimum  Space-Time  Processor 

4.1  Historical:  The  Displaced  Phase  Center  Antenna  (DPCA) 

4.1.1  The  DPCA  technique 

DPCA  (displaced  phase  centre  antenna)  is  a technique  which  compensates  physically  for  the  motion  of  the 
radar  platform  to  reduce  the  effect  of  motion  induced  Doppler  spread  of  clutter  returns.  Consider  two  antennas 
in  sidelooking  configuration  as  shown  in  Figure  22.  At  time  m = 1 they  assume  the  dashed  position,  at  time 
m = 2 the  solid  one.  As  can  be  noticed  the  first  antenna  at  in  = 1 assumes  the  position  of  the  second  sensor 
at  time  m = 2.  This  is  equivalent  to  having  one  antenna  fixed  in  space6  for  the  duration  of  one  pulse  interval. 
Clutter  suppression  is  done  by  subtracting  so  that  the  clutter  remainder  becomes  e = q2  — C21.  If  C12  = C21 
perfect  cancellation  is  obtained.  Notice  that  this  technique  was  implemented  in  RF  technology  long  time  before 
the  age  of  digital  signal  processing. 

4.1.2  A note  on  DPCA  and  STAP 

For  more  than  3 decades  the  DPCA  principle  (motion  compensation  by  spatial  coincidence  of  sensor  positions, 
s.  section  4.1.1  and  Figure  22)  has  been  considered  the  physical  background  of  space-time  clutter  rejection. 
More  recently,  numerical  investigations  have  shown  that  forward  looking  arrays  (which  do  not  have  the  DPCA 
property)  work  as  well  with  space-time  adaptive  processing. 

This  reveils  that  the  function  of  STAP  is  not  based  on  a DPCA  geometry.  The  space-time  filtering  is  just 
based  on  the  fact  that  airborne  clutter  echoes  are  signals  depending  on  the  two  variables  space  and  time,  and 
they  are  bandlimited  in  the  Doppler  as  well  as  in  the  azimuth  dimension.  If  such  signals  are  properly  sampled 
in  space  (sensor  displacement)  and  time  (PRF)  any  kind  of  filtering  can  basically  by  applied  without  aliasing 
losses.  The  property  of  slow  target  detection  is  based  on  the  special  shapes  of  such  clutter  spectra  (narrow 
ridge).  STAP  is  not  based  on  DPCA.  DPCA  is  merely  a special  case  of  STAP.  The  DPCA  property  plays  a role 
in  the  context  of  compensating  for  the  effects  of  system  bandwidth. 

4.2  The  LR-Test  for  2-D  Vector  Quantities 

The  principle  of  detecting  a signal  vector  s before  a noisy  background  given  by  q is  briefly  summarized.  Let 
us  define  the  following  complex  vector  quantities: 


( qi  \ 

( Si  \ 

( Xi  > 

qz 

S2 

X2 

q = 

; s = 

; 

; * = 

\ <1*  ) 

\ SN  / 

\ 

where  qm,sm  and  xm  are  the  spatial  subvectors  (signals  at  the  array  output)  at  the  m-th  pulse  repetition  interval. 
In  general  the  noise  vector  consists  of  a correlated  part  c (e.g.  jammer,  clutter)  and  an  uncorrelated  part  n (e.g. 
receiver  noise): 

q = c + n (18) 

The  signal  vector  s is  assumed  to  be  deterministic,  x is  the  actual  data  vector  which  may  be  noise  only 
(x  = q)  or  signal-plus-noise  (x  = q + s).  The  problem  of  extracting  s optimally  out  of  the  background  noise 
q is  solved  by  applying  the  well-known  linear  weighting 

wopt=7Q~1s  (19) 

6Due  to  the  factor  of  2 in  the  Doppler  term  2 mvpT  the  antenna  motion  during  one  PRI  is  only  half  the  antenna  spacing.  It  is 
important  that  phase  coincidence  of  clutter  echoes  occurs. 
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A block  diagram  of  the  optimum  processor  is  shown  in  Fig.  23.  The  spatial  dimension  is  given  by  the  N antenna 
elements  while  the  temporal  dimension  is  given  by  shift  registers  where  M subsequent  echoes  are  stored.  These 
space-time  data  are  multiplied  with  the  inverse  of  the  space-time  adaptive  clutter  covariance  matrix  for  clutter 
cancellation.  The  output  signal  are  then  fed  into  a space-time  weighting  network  whose  coefficients  form  a 
space-time  replica  (beamformer  and  Doppler  filter)  of  the  desired  signal. 

4.2.1  Performance  of  the  Optimum  Processor 

The  efficiency  of  any  linear  processor  w can  be  characterized  by  the  improvement  factor7  which  is  defined  as 
the  ratio  of  signal-to-noise  power  ratios  at  output  and  input,  respectively 


IF  = 


pin 

pin 


w'ss^w 

w'Qw 

s^g 

trtQj 


w*ss*w  • tr(Q) 
w’Qw  • s*s 


(20) 


Figs.  24  and  25  show  examples  for  the  improvement  factor  in  the  azimuth-Doppler  plane  for  sidelooking 
and  forward  looking  linear  arrays.  Along  the  clutter  trajectory  we  have  now  a clutter  notch. 


4.2.2  Comparison  with  1-dimensional  Methods 

The  question  is,  how  much  is  the  advantage  of  space-time  processing  versus  conventional  techniques.  Such  a 
comparison  has  been  made  in  Fig.  26.  The  optimum  space-time  processor  is  compared  with  a beamformer  cas- 
caded with  an  optimum  temporal  clutter  filter,  and  simple  Beamforming  plus  Doppler  filtering.  The  advantage 
of  space-time  processing  is  obvious. 


4.2.3  Range-Doppler  Matrix 

Plotting  clutter  power  or  IF  versus  Doppler  and  range  results  in  the  range-Doppler  Matrix.  Figures  27  and  28 
show  examples  for  sidelooking  and  forward  looking  . For  sidelooking  radar  (Figure  27)  the  clutter  trajectory 
is  a straight  vertical  line  which  corresponds  to  the  fact  that  the  clutter  Doppler  is  range  independent,  compare 
with  Figure  7.  For  forward  looking  radar  we  notice  a certain  dependency  of  the  clutter  Doppler  with  range, 
especially  at  short  range,  which  is  consistent  with  Figure  8. 

4.3  Optimum  Processor  and  Eigencanceller 

The  Eigencanceler  is  a zero  noise  approximation  of  the  optimum  processor.  It  is  given  by 

P = I — E(E*E)-1E*  (21) 

where  E is  the  matrix  of  eigenvectors  belonging  to  the  interference  component  in  Q.  Figure  29  shows  a 
comparison  of  the  optimum  processor  and  the  eigencanceller.  The  curves  are  almost  identical,  except  for  the 
clutter  notch.  Here  the  optimum  processor  suppresses  the  interference  down  to  the  noise  level  whereas  the 
eigencanceler  forms  an  exact  null.  The  eigencanceler  needs  less  training  data  for  adaptation  than  the  optimum 
processor  which  can  be  a significant  advantage  for  real-time  operation. 


7The  expression  improvement  factor  is  commonly  used  for  characterizing  temporal  (i.e.,  pulse-to-pulse)  filters  for  clutter  rejection. 
The  same  formula  may  be  used  lot  spatial  applications  (array  processing)  in  the  context  of  interference  or  jammer  suppression.  Then 
instead  of  improvement  factor  the  term  gain  is  used.  Since  our  main  objective  is  clutter  rejection  we  prefer  the  term  improvement  factor 
or  its  abbreviation  IF. 
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m — 1 

m = 2 


Figure  22:  Principle  of  a 2-pulse  DPCA  clutter  can- 
celler 


1 2 N 


Figure  24:  Improvement  factor  for  sidelooking  array 


Figure  25:  Improvement  factor  for  forward  looking  ar- 
ray 


Figure  23:  The  optimum  adaptive  space-time  proces- 


sor 
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Figure  26:  The  potential  of  space-time  adaptive  pro-  , . . . , . , 

b . Figure  28:  Range-Doppler  Matrix  (greytones  denote 

cessing:  o optimum  processing;  * beamformer  + IF/dB,  R = range/m,  FL,  <z>£  = 0°) 

Doppler  filter;  x beamformer  + adaptive  temporal  fil- 


-V7  „ . , Figure  29:  Comparison  of  optimum  and  orthogonal 

Figure  27:  Range-Doppler  Matnx  (greytones  denote  ...  & 


IF/dB,  R = range/m,  SL,  tpL  = 90°) 


projection  processing 
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